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Abstract 

The GRAPPA (GeneRalized Autocalibrating Partially Parallel Acquisi- 
tions) method of parallel MRI makes use of an autocalibration scan (ACS) 
to determine a set of synthesis coefficients to be used in the image recon- 
struction. For EPI time series the ACS data is usually acquired once prior 
to the time series. In this case the interleaved i?-shot EPI trajectory, where 
R is the GRAPPA reduction factor, offers advantages which we justify from 
a theoretical and experimental perspective. Unfortunately, interleaved R- 
shot ACS can be corrupted due to perturbations to the signal (such as 
direct and indirect motion effects) occurring between the shots, and these 
perturbations may lead to artifacts in GRAPPA-reconstructed images. Con- 
sequently we also present a method of acquiring interleaved ACS data in 
a manner which can reduce the effects of inter-shot signal perturbations. 
This method makes use of the phase correction data, conveniently a part of 
many standard EPI sequences, to assess the signal perturbations between 
the segments of i?-shot EPI ACS scans. The phase correction scans serve 
as navigator echoes, or more accurately a perturbation-sensitive signal, to 
which a root-mean-square deviation perturbation metric is applied for the 
determination of the best available complete ACS data set among multiple 
complete sets of ACS data acquired prior to the EPI time series. This best 
set (assumed to be that with the smallest valued perturbation metric) is 
used in the GRAPPA autocalibration algorithm, thereby permitting con- 
siderable improvement in both image quality and temporal signal-to-noise 
ratio of the subsequent EPI time series at the expense of a small increase 
in overall acquisition time. 
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1. Introduction 



1998: Pruessmann et aL, 1999; Griswold et al 



P arallel accelerated MRI (ISodickson and Manning . 119971 : Uakob et al. 

20021 ) increases image acqui- 



sition speed by using the spatial sensitivity of each receiver coil in an array 
of receiver coils, in addition to the spatial encoding provided by the applied 
linear magnetic field gradients used in conventional (non-accelerated) imag- 
ing. Relative to conventional imaging this additional spatial information 
allows one to reduce the number of acquired phase-encoded (PE) lines of 
data, while maintaining a desired digital resolution and field-of-view (FOV), 
thereby accelerating the data acquisition. This increase in the image acqui- 
sition speed is usually stated in terms of the reduction factor R which is 
defined, for fixed nominal digital image resolution and FOV, as the ratio of 
the number of acquired PE lines for a conventional scan to the number of 
acquired PE lines for a parallel accelerated imaging scan. 

For the sake of clarity we make some definitions before continuing with 
this introduction. In keeping with terminology initiated in Griswold et al. 
(IGriswold et all |2002j) we will refer to the data sampled during accelerated 
imaging as the reduced set of PE lines and we will refer to the PE lines which 
are absent from the parallel accelerated imaging data set, but present in a 
conventional data set, as the missing set of PE lines. In addition we will 
refer to the union of the reduced set and the missing set as the nominal set 
of PE lines - equivalent to the set of PE lines acquired during conventional 
imaging - since it is the sampling characteristics of this set which determine 
the nominal digital image resolution and FOV. The right-hand side of Figure 
[TJ and the lower part of Figure [2] illustrate these definitions as they pertain 
to an EPI time series. 



* Corresponding author 
** Principal corresponding author 

Email addresses: sheltraw@berkeley.edu (D. Shcltraw ), binglis@berkeley.edu 
(B. Inglis ) 



Complete ACS Data Set Time Series Data 



r 

N a 

L 



-N. 



1/ 



1 N s ' 



reduced 
(solid) 



missing 
(dashed) 



XN ( 



epi 



Figure 1: 1-shot ACS acquisition: This figure illustrates the ACS and time series It- 
space data associated with a single 2D image slice for the R = 3 case. The bold arrow 
indicates temporal order of acquisition and the X symbol denotes multiple acquisitions 
of the data contained within the vertical square brackets. The solid horizontal lines in 
k-space denote the acquired data. A complete set of ACS data consists of N a PE lines 
which are collected during 1-shot EPI acquisition. In the time series data the reduced set 
and the missing set of PE lines are respectively denoted by the solid and dashed lines. 
There are N s PE lines in the nominal set of PE lines (nominal matrix size: N s x N s ). 



The GRAPPA method of parallel imaging (IGriswold et all |2002| ) has 
been of great interest since, unli ke other parallel imaging methods such as 



SENSE (Sensitivity Encoding) (IPruessmann et all Il999l ). it does not re- 



quire explicit knowledge of the receive fields for each element of the receiver 
array. Instead, the GRAPPA method uses the data from the receiver array 
in an autocalibration procedure which estimates a set of synthesis coeffi- 
cients used to synthesize the missing set of PE lines from the reduced set 
of PE lines over the set of receiver coils. Central to the autocalibration 
procedure is the acquisition of autocalibration scan (ACS) data consisting 
of a subset of the nominal set of PE lines. This subset, comprised of PE 
lines near the center of k-space, will be referred to as the complete set of 
ACS data. This complete ACS data set is used to calculate the synthesis 
coefficients over the set of receiver array channels. The left-hand side of 
Figure [1] and the upper part Figure [2] illustrate these definitions as they 
pertain to N a PE lines of a complete ACS data set acquired prior to an EPI 
time series for 1-shot and i?-shot interleaved ACS, respectively. 
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Figure 2: R-shot segmented ACS acquisition: This figure illustrates the ACS and time 
series k-space data associated with a single 2D image slice for the R = 3 case. The bold 
arrows indicate temporal order of acquisition and the X symbol denotes multiple acqui- 
sitions of the data contained within the vertical square brackets. The solid horizontal 
lines in k-space denote the acquired data. A complete set of ACS data consists of N a 
PE lines which are collected during R shots of a segmented EPI acquisition. In the time 
series data the reduced set and the missing set of PE lines are respectively denoted by the 
solid and dashed lines. There are N s PE lines in the nominal set of PE lines (nominal 
matrix size: N s x N s ). 



For GRAPPA EPI time series the ACS data is usually collected onc e 



prior to the acquisition of the time series data (jSchmiedeskamp et al.L 120100 



Both 1-shot and f?-shot interleaved fc-space trajectories ( jSchmiedeskamp et al. 



20101 ) have been used to acquire a complete set of ACS data. Figures [T] and |2] 



depict the k-space trajectories for 1-shot and Z?-shot interleaved segmented 
acquisition respectively, of complete ACS data sets for the R = 3 case. The 
justification for the use of the J?-shot as opposed to the 1-shot trajectory has 
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not been clearly established in the literature, nor has a robust method been 
established for dealing with the potential consequences of inter-shot signal 
perturbations (such as motion) during acquisition of a complete i?-shot ACS 
data set. 

The intent of this paper is two-fold: (1) To justify the use of the i?-shot 
interleaved EPI trajectory instead of the 1-shot trajectory from a theoretical 
and experimental perspective when in the presence of significant main mag- 
netic field inhomogeneity, and (2) To demonstrate a method of acquiring 
the i?-shot interleaved EPI ACS data which minimizes potential artifacts 
due to signal perturbations occurring between the R shots. We will not 
attempt to answer the question of under what circumstances (main field 
strength, susceptibility differences, FOV choice, and reduction factor) the 
artifact from motion may be greater than that from trajectory incompati- 
bility. Our hope is that with refinement of the method presented here this 
trade-off will simply not be in question. 

Although in this paper we focus on the direct effects of head motion - 
because we can control it to some extent - we note that ACS data acquired 
in a segmented manner can potentially be contaminated by any pertur- 
bation to the data that occurs on a timescale smaller than that required 
to sample a complete ACS data set. Such perturbations may include, for 
example, motion effects such as B changes due to subject chest motion, 
magnetic susceptibility changes due to motion, perturbations to the time se- 
ries steady-state when through-slice-plane motion occurs and signal spiking 
due to electrostatic discharge between cables. 

We note that motion can perturb GRAPPA EPI time series data by two 
means: (1) Motion during the ACS data acquisition can lead to an inaccu- 
rate estimation of the synthesis coefficients, thereby leading to artifacts in all 
images of the time series along with a potential degradation of tSNR (tem- 
poral signal-to- noise ratio); and (2) Motion between the time of the ACS 
data acquisition and the time at which any particular image of the time 
series is acquired may make the synthesis coefficients inappropriate since 
the imaged object may have moved into regions where insufficient signal 
existed du ring the ACS data acquisition - an incomplete spatia l sampling. 
Cheng flChend . l2010h has investigated the use of TGRAPPA (IBreuer et al.l 



20051 ). a method not yet commercially available, to provide new ACS lines 



for each volume of an EPI time series. This has the benefit of providing 
motion-uncontaminated ACS throughout the time series, but comes at the 
expense of reducing tSNR in the absence of significant subject motion, com- 
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pared to the performance of GRAPPA. Thus, for fMRI applications, it is 
to be expected that GRAPPA would be preferable to TGRAPPA, provided 
GRAPPA can be made more robust to motion. 

One of the methods in this paper involves assess ing motion contami- 
natio n of ACS data by the use of a navigator echo fjEhman and Felmled . 



19891 ) and a perturbation metric. Other investigators have studied methods 



by which to assess and reduce motion contamination of MRI data albeit 
in applications other than GRAPPA EPI time series. For example, Kim 
and Hu have used navigator echoe s to limit the effect s of motion in fMRI 
studies using t he FLASH sequen ce (IHu and Kiml . Il994l ) and the interleaved 
EPI sequence ( iKim et all Il996l ). 2 D navigators have been u sed in conjunc- 



tion w ith read-out segmented EPI ( iHeidemann et al.l . |2010| ) ( iNguyen et al. 



19981 ) for diffusion imaging, the navigator being used to determine whether 
a particular diffusion weighted ima ge should be re-acquired due to motion 
contamination. Holdsworth et. al. f lHoldsworth et all 120091 ) have also used 



a k-space entropy metric to assess motion corrupti on in read-out seg mented 



EPI for diffusion- weighted imaging. Law et. al. fjLaw et al.l . 120081 ) used a 



sliding window approach to update the coil sensitivity maps for the TSENSE 



(adaptive SENSE incorporating temporal filtering) (IKellman et al.l . l200l[ ) 
method. Here we extend the application of navigator echoes to assessing 
motion between ACS segments, a direct extension of previous ideas, albeit 
with a new application. 

The method presented here makes use of multiple i?-shot interleaved 
ACS EPI data sets and phase correction echoes, which are already part 
of most commercial EPI sequences, to assess the motion between ACS EPI 
segments and produce a complete set of ACS EPI data that is minimally cor- 
rupted by motion. The phase correction echoes are therefore doing double- 
duty for they will be used in their usual capacity to eliminate Nyquist 
ghosting and they will be used as navigator echoes to assess motion. This 
complete set of ACS segments, assessed to be minimally contaminated by 
motion, is then used to estimate the GRAPPA synthesis coefficients and 
hence synthesize the missing PE lines for the entire time series. With this 
redundant ACS scheme, i?-shot accelerated EPI time series can potentially 
be reconstructed with less artifact and greater tSNR in the presence of mo- 
tion. This suggests improvements for fMRI applications which we evaluate 
with a simple measurement of tSNR. 
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2. Theory 



In the GRAPPA method flOriswold et all 120021 1 the missing set of PE 



lines, of a nominal set with inter-sample distance Ak y = 1/FOV in the PE 
direction, are synthesized according to the following equation: 

N c b 

S n (m + r, X) = 2j S ' n ( m + m ' Rl X ) W nn'm>r (1) 

n'=l m'=—b 

where R is the reduction factor, m = . . . — 2R, —R,0,R, 2R, . . . enumerates 
the reduced set of PE lines, r = 1, . . . , R — lisa PE offset from a PE line 
of the reduced set to a neighboring PE line of the missing set, N c is the 
number of coil elements in the receiver array, W nn > m ' r (T y ) are the synthesis 
coefficients and b fixes a finite number of local PE lines of the reduced set 
to be used in the synthesis of the missing PE lines. Equation [T] is used in 
conjunction with the sampled ACS lines to perform an autocalibration step 
by which one obtains the synthesis coefficients through a fitting algorithm. 
Once the coefficients W nn ' m > r are determined then Equation [1] is used to 
synthesize the missing set of PE lines from the reduced set of PE lines 
for each receiver coil. The images for each coil are then reconstructed by 
a simple FFT (Fast Fourier Transform) of the nominal s et of PE lines 



followed by a square-root of a sum-of-squares combination (IRoemer et al. 



19901 : lLarsson et all 120031 ) of the individual coil images to produce the final 
image. 

The number of ACS lines N a (see Figures [T] and [2]) is usually signifi- 
cantly smaller than the nominal size N s of the sampling matrix in the PE 
direction. For non-EPI imaging sequences the ACS data and the reduced 
data can usually be acquired together by adding a few additional k-space 
lines (for ACS) to the sequence. This method of collecting ACS data is usu- 
ally considered unacceptable for EPI time series because it would extend 
the echo train length significantly and require unacceptable compensatory 
adjustments to the image-space spatio-temporal sampling, and because it 
would exacerbate the EPI-inherent effect of geometric distortion due to main 
field inhomogeneity. Therefore, for GRAPPA EPI time series the ACS are 
usually acquired once prior to the time series and a single estimation of 
the synthesis coefficients is done followed by repeated calculations of the 
missing PE lines from the reduced data sets at each discrete time point n 
of the time series. 

There is some flexibility in the scheme used to acquire ACS for an EPI 
time series. For i?-fold acceleration it is possible to use 1-shot to i?-shot 



(interleaved) ACS trajectories, each with distinct advantages and caveats. 
For example, a rapid estimation of GRAPPA coefficients may be obtained 
with a 1-shot ACS (see Figure ED), using Ak matched to that of the n ominal 



set of PE lines fjPorter and Heidemannl . 120091 ; iHeidemann et all l2010l ) . This 



ACS data set should be relatively uncontaminated by sample motion during 
the ACS acquisition itself. However, the i?-shot trajectory (see Figure[2]) has 
the advantage of eliminating artifact that results when the signal dynamics 
over a 1-shot ACS sampling trajectory differ significantly from that over 
the reduced dat a set trajecto ry, as may occur in the presence of main field 



inhomogeneity (jChengi . 120101 ) or T 2 * signal relaxation. Appendix A shows, 
from a mathematical frames perspective, that the GRAPPA equations are 
expected to break down when significant main magnetic field inhomogeneity 
is present unless the ACS data is acquired in a matched i?-shot trajectory. 
This breakdown of the GRAPPA equations is expected to be exacerbated 
when using strong main fields since main field inhomogeneity increases with 
main field strength and since a large reduction factor is often used to control 
the resulting geometric distortions due to the field inhomogeneity. 

In the remainder of this paper we: (1) Demonstrate the effects of ACS 
trajectories unmatched with respect to field inhomogeneity perturbations 
and T 2 * signal relaxation; and (2) Propose and demonstrate a simple method 
by which residual aliasing due to motion between i?-shot ACS segments may 
be reduced. The phase correction lines (see Figure |3]) which precede each 
2D slice in commercial EPI sequences are used as navigator echoes to assess 
whether signal changes due to motion or other perturbations have occurred 
between the segments of a complete ACS data set. Using the data from 
these navigator echoes the following metric M s will be used to assess the 
perturbations between the R segments of the s th complete ACS data set: 



M, 



Hm=l 2n=l 2i=l 1 1 'Ssjm M S ; 



R-l 



s,j+l,m 



Ml 



(2) 



where S , S j m [ra], j, m and n are the FFT of a line of phase correction data, 
the segment index (j = 1 , . . . , R) , the coil element index and the sampling 
index in the frequency-encoding direction, respectively. We expect that 
complete ACS data sets with the smallest associated M s will produce EPI 
images with the least residual aliasing. 
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Figure 3: Sequence diagram for the acquisition of one complete set of ACS data. G x , G y 
and G z denote the frequency-encoding, phase-encoding and slicing gradients respectively. 
The gray-filled ADC regions identify the phase correction data lines which were also used 
as navigator echoes to detect motion between interleaved ACS segments. Note that the 
prephasing G y gradient is incremented for each ACS segment to achieve the interleaving. 

3. Methods 

A single phantom experiment was performed to demonstrate the arti- 
facts encountered when the 1-shot ACS trajectory is used in the presence 
of significant main field inhomogeneity. Two brain experiments were per- 
formed to show the efficacy of the proposed navigator-based method of 
acquiring ACS with reduced motion contamination. 

3.1. Data Acquisition 

All data were acquired on a Siemens TIM Trio 3T whole-body scan- 
ner with a 12-channel phased array head receiver coil. Brain images were 
obtained from volunteers in accordance with a protocol approved by the 
institutional Committee for the Protection of Human Subjects. 
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3.1.1. Phantom Experiment 

In the first experiment we acquired data to compare GRAPPA recon- 
structed EPI images resulting from f?-shot segmented interleaved ACS ver- 
sus 1-shot ACS. A structural phantom (The Phantom Lab, Salem, NY, 
USA) was used to assure that the image data was motion-free. An EPI 
sequence was modified in-house to either acquire ACS data for R = 3 
GRAPPA in a 1-shot acquisition or a segmented interleaved 3-shot acqui- 
sition. It should be noted that Siemens' commercial sequences use a 1-shot 
trajectory for R = 2 and an i?-shot segmented interleaved trajectory for 
R > 2 when acquiring the ACS data set. 

The relevant imaging parameters for this comparison were as follows: 
TR = 2450 ms, TE = 33 ms, nominal matrix size N s x N s = 96 x 96, FOV 
= 224 x 224 mm, number of ACS reference lines N a = 24, echo-spacing = 1.6 
ms and slice thickness = 3 mm. This experiment facilitated an evaluation of 
the potential for mismatched field inhomogeneity to affect residual aliasing. 
Note that an echo spacing roughly twice as long as is typical for fMRI was 
chosen to provide a clear example of the mismatch problem. We did not 
set out to evaluate the severity of the artifact for 1-shot f?-fold acceleration 
per se. 

3.1.2. Brain Experiments 

In a second experiment we demonstrate the proposed method of multiple 
ACS acquisition for eliminating residual aliasing due to motion during R- 
shot ACS. Multislice 2D brain images were obtained from volunteers using 
the following imaging parameters: GRAPPA with R = 2, echo-spacing = 
0.8 ms, nominal matrix size iV s x N s = 64 x 64, number of phase correction 
reference scans (which also serve as the navigator echoes) = 3, and TR = 
2000 ms. An EPI sequence was modified in-house to acquire N acs = 10 
complete R = 2 interleaved ACS data sets (see Figure [2] for a depiction 
of the analogous R = 3 case) instead of the usual single complete set. 
Each interleaved ACS EPI segment acquired 12 PE lines for a total of 24 
ACS lines (N a = 24) in a complete ACS data set. Following the ACS 
acquisition a single volume of 2D multislice image data (the reduced set) 
was acquired. For fMRI applications a time series of reduced data volumes 
would be acquired, but in this experiment we were only interested in the 
effects of motion during the ACS and therefore a single reduced k-space 
volume was sufficient. 

In order to investigate the various effects of subject motion during in- 
terleaved ACS, four separate trials were acquired in the second experiment. 
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During the acquisition of the first three trials the subject nodded their head 
(approximately 3 degrees) in randomly distributed 4-second intervals, be- 
tween which the subject tried to remain motionless. During the acquisition 
of the fourth trial the subject tried to remain motionless throughout the ac- 
quisition of the 10 complete ACS data sets. The fourth trial was considered 
to be the target, or "best case", data. 

In a third experiment we sought to assess the effect of motion contami- 
nated ACS on the tSNR of a time series EPI acquisition, as would be used 
for fMRI. Once, again the subject was instructed to nod their head during 
ACS acquisition after which they were to remain motionless throughout an 
EPI time series. The relevant imaging parameters were R = 2, number of 
time series volumes N epi = 100, nominal matrix size N s x N s — 64 x 64, 
number of ACS reference lines N a = 24, number of repetitions of segmented 
interleaved ACS complete data set N acs = 10 (as depicted in Figure EJ), 
number of phase correction reference scans (which also serve as the navi- 
gator echoes) = 3, TR = 2000, TE = 28 ms, echo-spacing = 0.52 ms and 
FOV = 224 mm. 

3.2. Image Reconstruction 

Image reconstruction was performed with either commercial Siemens 
GRAPPA algorithms and code or with in-house GRAPPA reconstruction 
code which implemented "method 1" (a PE-only GRAPPA interpolation 
with a sma l l num ber of blocks and no sliding window averaging) given in 



Brau et all (120081 ). 



3.2.1. Phantom Experiment 

For the first experiment image reconstruction was done on the scan- 
ner with standard Siemens reconstruction code for both the product and 
modified sequence. 

3.2.2. Brain Experiments 

For the second and third experiments the image reconstruction was done 
with in-house GRAPPA code. The EPI sequence acquires three phase cor- 
rection lines prior to each 2D slice of image data which we used as navigator 
echoes. Using Eq. [2] we calculate M s associated with each phase correction 
line and then we simply averaged the three M s values to yield the final M s 
value for each of the complete ACS data sets of a given trial. 

The four trials of the second experiment were processed in the manner 
of four "ACS time series" acquisitions, with each set of the ten interleaved 
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ACS being combined with the single reduced EPI data set following each 
looped ACS. The 20 separate ACS segments from ten interleaved R = 2 
ACS acquisitions can be grouped into 19 complete ACS data sets using 
adjacent pairs only. We do not expect the order of the interleaved ACS 
acquisition to matter, but we chose to restrict the metric M s to ACS pairs 
acquired as nearest neighbor pairs in order to minimize motion and scanner 
drift effects. From each complete ACS data set the synthesis coefficients 
were estimated. The missing set of PE lines was then synthesized using the 
reduced data set of PE lines and the synthesis coefficients corresponding to 
each of the N acs = 19 complete ACS sets. This was done for each of the 
twelve coil elements, after which the data was Fourier transformed to image- 
space and the images from each coil were combined in a final sum-of-squares 
magnitude image. 

Data from the third experiment was processed in two ways. As in the 
second experiment, M s was evaluated for 19 nearest neighbor pairs of com- 
plete ACS data sets from the 10 acquired ACS pairs. From these, the ACS 
pairs with the highest and lowest M s were determined and used in the re- 
construction of the same N epi = 100 volume time series EPI data, from 
which tSNR maps were produced. TSNR images were obtained by calcu- 
lating the ratio of the temporal mean to the temporal standard deviation 
over the EPI time series for each pixel in image space. 



4. Results 

4-1. Phantom Experiment 

The bottom of row of FigureHJshows, for the case of R = 3, an example of 
the residual aliasing artifact that may result when 1-shot ACS data is used to 
autocalibrate the synthesis coefficients. Such artifact is clearly undesirable 
and is much reduced by using 3-shot segmented interleaved ACS instead, 
as shown in the top row of Figure HI We did not seek to quantify the 
extent of residual aliasing because the amount depends on R, the main field 
inhomogeneity (which itself may depend upon main field strength and the 
heterogeneity of the magnetic susceptibility of the imaged object), echo- 
spacing time and the FOV. The effects of field inhomogeneity upon the 



GRAPPA equations are considered in Appendix A 



12 



j^to^. jt^^b\ ^wl Ifct 



^WHft^, j^tBpv yw-Bry 



Figure 4: Top row: Five successive slices of R = 3 GRAPPA images reconstructed using 
a 3-shot interleaved ACS data acquisition. Bottom row: The same set of image slices 
but using a single-shot ACS data acquisition. Note the residual aliasing (white arrows) 
which extends throughout the images of the bottom row. Exterior to the phantom the 
residual aliasing appears as misplaced image intensity, while interior to the phantom the 
residual aliasing appears as both increased and decreased local image intensity. 



4-2. Brain Experiments 

Figures [5] and [6] show the results for a GRAPPA-EPI R = 2 acquisition 
acquired in the second experiment. Figure [5] shows plots of the perturbation 
metric M s associated with each of the 19 complete ACS data sets for the 
four motion trials. Although the perturbation metric in the first three trials 
attains values nearly an order of magnitude greater than those of the fourth 
(no intentional motion) trial, it is clear that low M s sets of ACS data are 
available in all cases. 

Figure [6] shows the EPI images associated with successive complete ACS 
data sets. For brevity only the 10 images associated with the even valued 
complete ACS data sets are shown, rather than all 19. Image rows 1 through 
4 correspond to the four different trials during which varying amounts of 
motion were introduced. Image row 4 was acquired as the subject tried to 
remain motionless throughout ACS acquisition. For each of the runs, the 
dashed boxes enclose the GRAPPA reconstructed images for which M s was 
highest and the solid boxes enclose the images for which M s was lowest. In 
each of the enclosed images the value of M s is given in the upper left corner. 
Figures [5] and [6] show a clear correlation between the perturbation metric 
M s and the visually apparent residual aliasing. 

Figure [7] shows two mosaics of tSNR images for the same 100 volume 
EPI time series data acquired in the third experiment. The top mosaic 
is associated with the largest valued perturbation metric during the ACS 
while the bottom mosaic is associated with the smallest valued perturbation 
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Figure 5: Plots of the perturbation metric M s as a function of the n th complete ACS 
data set. Plots A through D correspond respectively to rows 1 through 4 of Figure |6] 
Note that the scale of the ordinate varies between the four plots. 

metric. Both images were constructed using the same gray-value scale. 
Visual inspection of these images shows significant improvement in tSNR of 
the EPI time series when the ACS data associated with the smallest valued 
perturbation metric is used in the calculation of the GRAPPA synthesis 
coefficients. For reference note that for the 5x5 pixel region shown in Figure 
[7] the ratio of the mean tSNR for the largest valued metric case to the mean 
tSNR for smallest valued metric case is tS N R sma u est /tS N Ri argest = 1.54, 
the improvement in tSNR being due to the proper localization of the signal 
in image space. 
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Figure 6: Images reconstructed from ACS data sets during which subject motion was 
introduced. Image rows 1 through 4 correspond to four different trials. During the ac- 
quisition of trials 1-3 the subject was instructed to move during ACS acquisition. During 
trial 4 the subject was instructed to remain motionless during the ACS acquisition. For 
brevity only the 10 images associated with the even valued complete ACS data sets are 
shown rather than all 19. For each trial the image corresponding to the highest or lowest 
value of the ACS perturbation metric M s is indicated by a dashed or solid bounding box 
respectively. The value of M a is given within the bounding box. 



5. Discussion 

In this work we demonstrated the effect of field inhomgeneity for one 
example situation to highlight the severity of artifacts that can result from 1- 
shot ACS. (The mismatch phenomenon is treated in depth in Appendix A| ) 



With increasing R, field inhomogeneity or EPI echo-spacing time the effects 
of main field inhomogeneity on 1-shot ACS data acquisition are expected to 
become increasingly pronounced, suggesting a growing need for interleaved 
ACS data acquisition when acquiring ACS data for GRAPPA-accelerated 
EPI at 3T and above. 

It is important to not substitute one source of error for another when- 
ever possible. Therefore we have developed and demonstrated in this paper 
a means of acquiring segmented interleaved ACS data which greatly dimin- 
ishes ACS motion contamination. Acquiring redundant ACS prior to an 
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Figure 7: 2D multislicc mosaics of tSNR images for the EPI time series. Top: Time series 
reconstructed by using the complete ACS data set with the largest valued perturbation 
metric. Bottom: Time series reconstructed by using the complete ACS data set with 
the smallest valued perturbation metric. Note that the same 100 volumes of accelerated 
(R = 2) time series EPI data were used here, and only the ACS used to compute the 
synthesis coefficients differs between the two sets of images. For reference, the ratio of 
the mean tSNR, over the 5x5 pixel region outlined in white, for smallest and largest 
valued M s is tSNR smaUest /tSNRi arges t = 1.54. 



accelerated EPI time series allows the use of navigator echoes, conveniently 
already part of a product EPI pulse sequence, as a means of assessing mo- 
tion corruption of GRAPPA ACS data. We have demonstrated that it is 
possible to retrieve the lowest motion-contaminated ACS data from a series 
of ACS data sets. Moreover we have demonstrated that motion during ACS 
has important consequences with respect to image artifact and time series 
tSNR. It is expected that this simple modification will greatly enhance the 
robustness of GRAPPA-accelerated EPI for fMRI, where subject motion 
during the initial part of the acquisition - the interleaved ACS - can render 
an entire time series worthless in a worst-case scenario, or with compromised 
tSNR as demonstrated in Figure UJ 

In this work we acquired a fixed number of complete ACS data sets but 
the method could be further extended by looping on the ACS acquisition 
until the metric M s meets a prior criterion, thereby assuring that acceptable 
GRAPPA coefficients will be obtained in the presence of protracted motion 
during ACS acquisition. Figures |5] and [6] show that it should be possible to 
establish such a criterion for the metric M s , although it would necessitate an 
automated way to advance the pulse program to the time series acquisition 
phase. 

One limitation of the present method is a potential directional depen- 
dence of the navigator echoes. Since the navigators are projections of a 2D 
image slice (multiplied by a receive coil field) onto the frequency encoding 
axis, this method should be relatively insensitive to motion that is in the 
PE direction only. The nonuniformity of the receive coil field will introduce 
some sensitivity to motion that is in the PE direction only but if enhanced 
sensitivity to this type of motion is needed then a second navigator, spatially 
orthogonal to the first, could be used. 

The work reported in this paper stands in con trast to the T GRAPPA 
EPI time series method investigated by Cheng ( Cheng . 2010l ). In the 



TGRAPPA EPI times series method the effects of motion upon autocal- 
ibration are addressed by updating the synthesis coefficients at each point 
in the time series. In our work we take the perspective that a single cal- 
culation of the synthesis coeffcients will generate the greatest stability in 
the synthesis of the missing PE lines. We find support for this perspective 
in Cheng's TGRAPPA work where he noted that in the absence of signif- 
icant motion the use of frequently updated synthesis coefficients leads to 
a reduction of tSNR compared to the usual single calculation of s ynthesis 
coefficients. Indeed, in separate work, Cheng has shown recently (jChengi . 
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20121 ) that the temporal noise properties of GRAPPA EPI time series can 
be reduced via a single set of ACS for multiple time series acquisistions, 
rather than one ACS per acquisition. This increases still further the impor- 
tance to acquire a motion- free ACS. Thus, the method presented here shows 
that it is possible to reduce the effects of motion during ACS acquisition 
while obtaining the desirable properties of the i?-shot interleaved ACS, i.e. 
matched distortion characteristics for the ACS and the accelerated time 
series data. This should be of special interest in fMRI where motion is 
especially problematic. 

From our perspective any motion during the time series and following 
the ACS leads to a spatial undersampling of the receive coil fields - the head 
may have moved into a region of the array from which no signal (or signal 
of low SNR) originated during the ACS. We see this spatial undersampling 
as a separate issue to be addressed in future work. It is a more involved 
problem than the simple method presented here that permits the acquisition 
of ACS uncontaminated by motion. 

Appendix A. GRAPPA and Field Inhomogeneity 

We assume that there are N c receive coils in the receive array. For the 
reduced data set the signal from the nth receiver coil, after Fourier transform 
in the f req u ency-enco ding x-dir ection and phase correction (Nyquist ghost 



correction) (IReeder et all Il999l ) , may be written as 



where 



(x, mR) = (p, g nm ) = J p(x, y)g* m {x, y)dy (A.l) 



g nm (x, y) = g* n (x, y ) e »*«^ e -^.«nv) = E m g* n {x, y). (A.2) 



and where the definition of E should be obvious, g n (x, y) is the receiver coil 
field (coil sensitivity), R (an integer) is the reduction factor, <fi(x,y) is the 
local field inhomogeneity offset of the main field, r e is the echo-spacing time 
and meZ. Note that p(x, y) is assumed to be an image that is compactly 
supported on |y|,|x| < 1/Afc and that the PE direction, y, will be the 
direction of acceleration. Also note that s n , g n and </> will depend upon z 
as well but we have omitted explicitly writing this dependence. Note that 
although (p, g nm ) is a function of x we will suppress the x-dependence for 
the sake of an economy of symbols. 
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When the ACS data are acquired in R interleaved segments the signal 
in the n th receiver coil will be 



s%(x, mR + p) = (p, g* m ) = J p(x, y)g^(x, y)dy (A.3) 
where p = 0, . . . , R — 1 and 

ffijx, y) = 9*(X, y yMmR +P )Aky e -imr^,y) = ETE^g^X, y). (A.4) 

The superscript R denotes the R-shot trajectory for the complete ACS data 

set. When the complete ACS data set is acquired in one segment (1-shot 
trajectory) then the signal may be written as 

40, mR + p) = (p, g l nm ) = p(x, y)gl l * m (x, y)dy (A.5) 



where 

gUx,y) = g * {Xj yyMrnR +P )Aky e -i(rnR +P )r^ = E™ E{ o g* n {x , V ) . (A.6) 

and where the definition of the modulation operators Er, Er E\ and E\ 
should be obvious. The superscript 1 denotes the 1-shot trajectory for the 
complete ACS data set. 

For a suitable g n (x, y) and R the set {qw l, where n = 0, . . . , N c — 1 and 
meZ, will form a frame ( Christensenl . 120031 ) . Since {g nm } (reduction factor 



R > 1) is assumed to be a frame we can write 

oo JV C 



p= 2 Yj^ p ' gm '^ s lgn ' m ' ( A - 7 ) 

m'=— oo n'=l 



where S is the frame operator for the frame \q nm } and where \S 1 g nm } is 



the canonical dual frame to the frame {g nm } ( Christensen . 20031 ). Note that 



the dual frame depends upon x although we will usually not indicate this 
explicitly. Note also that the canonical dual frame is only one of the possible 
and non-unique dual frames for which Equation f]A.7|) may be written. The 
canonical dual frame is a special dual frame in that it has the property that 
the (p, S~ l g m ^ n iy have a minimal / 2 -norm among all possible dual frames. 
If the inverse frame operator could be found then equation f lA.7j) would 



form the fundamental means of obtaining p. In practice obtaining the in- 
verse frame operator is mathematically difficult to obtain and this is at 
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least partly due to the lack of an explicit knowledge of the coil sensitivities 
g n (x,y) (or the field inhomogeneity which may perturb the frame). Au- 
tocalibration methods provide an alternative to obtaining estimates of the 
dual frame and instead focus on determining the coefficients that trans- 
form one frame to another. We now take a closer look at the GRAPPA 
autocalibration method and how main-field inhomogeneity may affect the 
method. 

Performing an analysis of p given by equation flA.Tj) with respect to the 
frame {<3^ m } we obtain 



oo N c 

(P,9nm>= 2 Yl(Pi9 n ' m 'X S ~ 1 9n>m>,giL)- ( A ' 8 ) 
m'=— oc n'=l 

which we may write as 

oo N c 

(P,9L)= 2 'Z<p i g n 'fnfXS- 1 E mf g^ES^ to g n y. (A.9) 

m'=— oo n'=l 

Since the frame operator S, and there f ore it s inverse, commutes with the 
modulation operator E m (jChristensenl . 120031 ) then equation (IA.80 may be 
written as 

oo N c 

(P,9L)= 2 Y,^9n' m '){E m '- m S- l g:„El g n ). (A.10) 

m' =— oo n'=X 

or according to equations (lA.ip and (1A.3j) 

oo N c 

s^(x,mR + p)= ^ s n '(x,m' R) c(m' ~ m,n' ,n,p, x) (A. 11) 

m'=— oo n'=l 

where we have defined the synthesis coefficients 

c{m'-m,n\n,p,x) = (E^S^g*,, E p m g n ) (A.12) 

Notice that the c depend upon the difference m — m! only. This particular 
dependence upon m and m! allows us, through a change of indices and a 
truncation of the sum over m! (which must occur in a any practical setting), 
to write 

b 2 N c 

Sn(x,mR + p) = ^ s n /(x,mR + m' R) c(m' ,n' ,n,p,x) (A.13) 

m! =b\ n'=l 

20 



which is equivalent to the form of the GRAPPA equations given by Equation 
([T]). The quantity B = b 2 — b\ is often referred to as the block size. 

When the GRAPPA equations and the complete ACS data set are used 
to solve for the coefficients c(m' ,n' ,n,p, x) it is usually the case that the 
number of unknown coefficients c(m', n', n,p, x) exceeds the number of equa- 
tions. When the number of ACS lines N a is an interger multiple of the re- 
duction factor R then there will be N c N a (R— 1)/R equations and BN c N a /R 
unknown coefficients c(m', n', n,p, x) at each point x. For a typical example 
of N a = 24, N c = 12 and R = 2 the number of equations would be 144 and 
the number of unknowns is 1445. This presents a problem with respect 
to meeting necessary (but not sufficient) requirements for a unique solution 
set of the c(m',n',n,p,x) - number of equations must be greater than or 
equal to the number of unknowns. To meet the necessary requirement it 
is usually assumed that the c(m',n',n,p,x) vary little with respect to x in 
the neighborhood of any x in the FOV but that the signal does vary signif- 
icantly with x. This allows one to use the ACS data at neighboring x (at 
least (d 2 — d 1 )/(R — 1) are needed) to locally determine the c(m',n',n,p,x). 
Such an assumption implies that g n , S~ l and <fi vary little about any x in 
the FOV. 

When the 1-shot trajectory is used to sample the complete ACS data 
set then, following steps similar to those given above, we can write 

oo N c 

s l n {x,mR + p)= J] J^s n ,(x,m , RXE m 'E; m S- 1 gl,,E% g n ). (AAA) 

m'=—co n' = l 



Since E E\ it is not possible to cast Equation ( 1A.14|) into the form of 



Equation ( 1A.11)> because the indices m and m' do not appear as a difference 
only. Therefore, except when (f)(x, y) = 0, it is necessary that the ACS 
data be collected over an R-shot trajectory in order for GRAPPA equations 
to strictly apply. As R or <p(x, y) increases in magnitude the departure 
from the GRAPPA equations will increase. We can therefore expect the 
artifacts resulting from a 1-shot ACS calibration scan to increase with main 
field strength where susceptibility effects will be exacerbated and increasing 
reduction factors may be applied to combat the resulting distortions due to 
field inhomogeneity. 
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